Scaling of random spreading in small world networks 
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In this study we have carried out computer simulations of random walks on Watts-Strogatz-type 
small world networks and measured the mean number of visited sites and the return probabilities. 
These quantities were found to obey scaling behavior with intuitively reasoned exponents as long as 
the probability p of having a long range bond was sufficiently low. 
PACS numbers: 05.40.-a, 05.50.+q, 87.18. Sn 



Recent interest in some disordered graphs was stimulated 
by a series of discoveries related to the so called small 
world properties of certain networks including social nets 
document retrieval in the www §, the internet g, 
scientific cooperation relations Q etc. (for reviews see 
1^). Here 'small world' is used in the sense, that arbi- 
trarily selected pairs of nodes can be reached from each 
other within few steps on the average, in spite of the 
relatively low number of links present in the system. 

Many of the earlier network studies have dealt with the 
Watts-Strogatz small world network (WSSWN) model, 
which we will also study. We define it as a one- 
dimensional ring consisting of N nodes with fc-neighbor 
interaction where, in addition, N ■ p new links are in- 
troduced between arbitrarily chosen, not yet connected 
nodes. Though the small world property is closely re- 
lated to the dynamic process of information spreading in 
the network, relatively few papers have been devoted to 
this aspect [^-^. 

In the recent paper Jasch and Blumen |^ studied the tar- 
get problem on a WSSWN model, i.e., the decay of the 
number of randomly distributed target sites on the graph 
where annihilating random walkers move around. There 
it turned out that this decay process is closely related to 
the mean number of distinct visited sites S{n), where n is 
the number of steps done by the walkers. In addition they 
observed the expected crossover in S{n) from the ^/n be- 
havior that is characteristic for the one-dimensional case 
to S(n) oc n describing the high dimensional or random 
graph situation. S{n) as a function of p and n has a 
scaling form: 
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where / is a universal scaling function with the following 
properties: 
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const for x ^ 1 
^/x for X ^ 1 
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of the average distance between nodes having long range 
links for which the walker needs (x steps to 
sweep through. Thus the argument of the scaling func- 
tion / in Eq. (1) should be n/n^. In spite of this strong 
argument Jasch and Blumen found in their numerical 
simulations of WSSWN's a value a = 1.85. In the simu- 
lations they had chosen N = 50000 by taking an average 
over 500 random walkers for each of the 100 WSSWN and 
they varied p in the interval 0.01 < p < 0.1. The main 
aim of the present note is to resolve the discrepancy in 
the value of a. 

It is first noted that Eq. (1), as is usual in scaling theory, 
is valid only asymptotically and in this case the scaling 
limit is A'^ ^ oo, n oo and p 0. The scaling regime 
can be estimated from the variation of the mean vertex 
distance ^ as a function of p 0: It turns out that £k/N 
is a scaling function with the argument x = p ■ k ■ N and 
of sigmoidal shape. This curve suggests that one can not 
expect good scaling for the above mentioned crossover, if 
p ■ k ■ N ^ 100. Therefore, it seems likely that in 1^ the 
investigated values of p were not small enough to assure 
the proper scaling behavior (p,„in ■ k ■ N — 1000 in |^). 

For this reason we have carried out our simulations with 
considerably smaller values of p. In order to do so, we had 
to increase the system size as well. Hence we have studied 
WSSWN's with k = 2, N ^ 10^ and for each value of 
p = IQ-'^, 10-3■^ 10-3, 10-2-5 we did 100 realizations and 
had 100 random walkers per reahzation {pmm-k-N = 20). 
In these simulations we measured the average number of 
distinct visited sites S{n) as functions of n and p, as 
depicted in Fig 1. In this plot it is seen that for the two 
largest values oip saturation of S{n) has set in. 



Intuitively it is expected that a = 2 since in the system 
there exists a basic length scale ^ oc 1/p, characteristic 
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FIG. 1. Raw data for the average number of distinct sites 
visited S as a function of time steps n and the probabihty 
p = 10"*, 10"^'^, 10"^, 10"^ '^ increasing from bottom to top. 
For large p the saturation due to the finite size A'' = 10^ of 
the systems is visible. 
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FIG. 2. Scaling plot of the data of Fig. 1 with a = 2. The 
inset shows the scaling plot with the exponent of [9]. 

Fig. 2 shows a scaling plot of the results on S{n) where 
S{n)/s/n is plotted vs. np". The scaling was found to 
be optimal with the ehoise of a = 2. For comparison, 
we have also shown the same plot with a = 1.85, the 
value found in . Our results clearly support the simple 
scaling picture discussed above. 

We have also checked the scaling behavior for another 
quantity, namely the return probability Poo('^)- This is 
known to decay as for the p = case while an 

exponential decay is expected for large p. A scaling form 
similar to Eq. (1) should be valid for Pqo- 

Poo = n-^/^ip{x) (3) 
where ip{x) is a rapidly decaying scaling function with 



the limit ip = const, for a; <C 1. However, the argument 
of if should be the same as in eq. (1-2), namely x = p"n. 
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FIG. 3. Raw data for the return probability Poo. The p 
values are the same as for Fig. 1, now increasing from top to 
bottom. The whole time interval was binned by 100 bins of 
equal sizes on the logarithmic scale. 




FIG. 4. Scaling plot of the data of Fig. 3 using a — 2. 
For minimizing the finite size effects the asymptotic value 
l/N — 10~* was subtracted from Poo. For comparison, the 
inset shows the scaling plot with the a of [9]. 

Here in the interest of even higher accuracy we have used 
10 times more runs for the averages. In order to minimize 
the effect due to the finite size of the samples, we have 
subtracted 1/A'^, i.e., the n —^ oo limit, from the measured 
values. Fig. 3 shows the data for the return probability 
Poo and Fig. 4 the scaling plots. Again, we see that the 
scaling with the intuitively expected a = 2 is superior to 
the one obtained by the value of 

In conclusion we have shown that for sufficiently small 
probabilities of long range links the proper scaling vari- 
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able for the average number of distinct sites visited and 
also for the return probability is np^ , i.e., the natural ex- 
ponent a — 2 holds for the Watts-Strogatz small world 
network. 
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